MacFormat 1994 August/August CD.bin/Shareware/Education/numericalmethods Folder/chap_5/lsline.m

Text File | 1994-06-05 | 373 b | 14 lines | [MATF/MATL]

function [a,b] = lsline(X,Y)
% [a,b] = lsline(X,Y)
% Construct the least squares line.
% X is an 1xn abscissa vector, input.
% Y is an 1xn ordinate vector, input.
% a is a coefficient; y = ax+b, output.
% b is a coefficient; y = ax+b, output.
xmean = mean(X);
ymean = mean(Y);
sumx2 = (X-xmean)*(X-xmean)';
sumxy = (Y-ymean)*(X-xmean)';
a = sumxy/sumx2;
b = ymean-a*xmean;